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Abstract. Symmetric nuclear matter is studied within the conserving, self-consistent T-matrix approxima- 
tion. This approach involves off-shell propagation of nucleons in the ladder diagrams. The binding energy 
receives contributions from the background part of the spectral function, away form the quasiparticle peak. 
The Fermi energy at the saturation point fulfills the Hugenholz-Van Hove relation. In comparison to the 
Brueckner-Hartree-Fock approach, the binding energy is reduced and the equation of state is harder. 
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The calculation of nuclear matter properties from the 
basic nucleon-nucleon interaction has been extensively stud- 
ied using Brueckner type resummation of ladder diagrams. 
This resummation allows to rewrite the ground-state en- 
ergy of nuclear matter using as an effective interaction, 
the G matrix, which takes care of the short range re- 
pulsive core in the nucleon-nucleon interaction [n] . Calcu- 
lations using realistic interactions lead to results, which 
lie along a line (the Coester line) shifted with respect 
to the phenomenological saturation point (po — .16fm~ , 
E/N ~ — 16MeV). The remaining discrepancy can be at- 
tributed to relativistic effects and three-body forces con- 
tributions 

The results on the binding energy depend on the single 
particle energies used in the kernel of the Bethe-Goldstone 
equation The so called standard choice uses a self- 
consistent auxiliary potential defined by the G-matrix be- 
low the Fermi energy and the free dispersion relation above 
kp- Another choice is to use the self-consistent poten- 
tial also above the Fermi momentum which gives the so 
called continuous choice for the single-particle energies 
in the Bethe-Goldstone equation. In Brueckner-Hartree- 
Fock (BHF) calculations the hole line expansion, irrespec- 
tive of the choice of the auxiliary potential, is believed to 
converge to values close to the BHF with the continuous 
choice for single-particle energies Q. 

Recently self-consistent approaches based on the in 
medium T-matrix approximation for nuclear matter have 
been studied In this way a spectral function 

for nucleons in nuclear matter including two-particle cor- 
relations is obtained. The ladder diagrams involved in the 
calculation of the in medium T matrix include also hole- 
hole propagation. The T-matrix approximation takes into 
account some of the higher order hole line contributions as 



compared to the G-matrix approach. It would be instruc- 
tive to study the saturation properties of nuclear matter 
for the self-consistent T-matrix approximation with real- 
istic interactions. 

The T-matrix approach is a ^-derivable approxima- 
tion The self-energy is constructed as a functional 
derivative of a set of two-particle irreducible diagrams. 
This assures the fulfillment of thermodynamical relations 
for the quantities obtained The most famous such a 
relation is the equality of the Fermi energy and binding 
energy at the saturation point 



Ef = E/N 



(1) 



The realization of the above relation is very important 
since it would give confidence to the single particle prop- 
erties obtained in the calculations. In Ref. |^ we studied 
the self-consistent T-matrix approximation with a simple 
interaction confirming to a very good accuracy the fulfill- 
ment of thermodynamical relations by the numerical solu- 
tions. In BHF calculations the Hugenholz-Van Hove rela- 
tion is badly violated. This discrepancy can be reduced by 
invoking rearrangement terms for the Fermi energy [ p^ , [l4| , 
[Tst . By construction, the single-particle energies obtained 
in the T-matrix approximation come out consistently with 
thermodynamical observables. Thus we expect that single- 
particle energies, scattering width or spectral functions di- 
rectly obtained from the self-consistent T-matrix approx- 
imation are meaningful |p^ . 

For attractive interactions cold nuclear matter forms a 
superfiuid. Calculations using dressed propagators in the 

ferfluid phase show a strong reduction of the gap |^,p^. 
We expect that around the saturation point the super- 
fluidity is very weak This means that the correction 
from the superfluid correlation energy to the binding en- 
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ergy is small. We restrict ourselves to normal nuclear mat- 
ter for all densities studied here. It allows us to compare 
with BHF calculations which are performed exclusively in 
the normal phase of nuclear matter. 

The results here presented are obtained using a sep- 
arable parameterization of the Paris potential p8| for S, 
P, D and F partial waves, for symmetric nuclear mat- 
ter. Wc use rank 3 and rank 4 parameterization for the 
and ^Si Di partial waves. In the ^Pq partial wave 
we use Mongan I interaction, in order to avoid unphysical 
resonances far off-shell. In the numerical iteration the full 
spectral function is discretized. For momenta close to the 
Fermi momentum the spectral function is separated into 
a background part and a quasiparticle peak approximated 
by a delta function. The numerical treatment of the en- 
ergy integrations for the spectral functions is done using 
convolution algorithms [p^|. 

The T-matrix approximation resumes ladder diagrams 
with dressed particle-particle and hole-hole propagators 



<p|T(P,f2)|p >^V{p,p) 

duji f dw2 f (Pq 
2^ 



+ 



(1 - /K) - 



27r J (27r)3 ' Q 
Aipi,u;i)A{p2,UJ2)<q\TiP,n)\p > 



(2) 



where pi,2 = P/2±q and /(w) is the Fermi distribution. 
The imaginary part of the corresponding retarded self- 
energy can be obtained closing a pair of external vertices 
in the T-matrix with a fermion propagator 



( 



duJi 

< (p-k)/2|Imr(p4 



-LUl)\{p-k)/2>A 



where 



{uj - p2/2m - ReE{p, w)) + lmS{p, uj^ 



(3) 



(4) 



is the self-consistent spectral function of the nucleon and 
b{uj) is the Bose distribution. The real part of the self- 
energy is related to ImS by a dispersion relation 



ReE{p,u)^SHFip)+r 



(5) 



with Shp{p) the Hartree-Fock self-energy. Eqs. |2[ ^, || 
and ^ are to be solved iteratively and at each iteration the 
chemical potential fi is adjusted to fulfill the condition on 
the density p 



d(jj 
2^ 



d^p 
(2^ 



Aip,u)f{uj) ^ p 



The spectral functions obtained in the self-consistent 
solution consist of a quasiparticle peak and a broad back- 
ground (Fig. |l|). As function of momentum the position 
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Fig. 1. The spectral function A{p,uj) as function of energy 
for p — 0, 255 and 340MeV (solid, dashed, and dotted lines 
respectively) . 



of the peak in the spectral function follows approximately 
the quasiparticle dispersion relation 



h ReS{p, ujp) . 

m 



(7) 



The quasiparticle peak is very sharp for momenta close to 
the Fermi momentum. It is a manifestation of the quasi- 
particle nature of excitations close to the Fermi surface. 
Indeed we find that the single-particle width — 2Imi7(p, ujp) 
is proportional to {p—pr)^ close to the Fermi momentum. 
The background of the spectral functions extend far from 
the quasiparticle peak. The part of the spectral function 
below the Fermi energy leads to nonzero occupancy for 
momenta above pp and gives a large, negative contribu- 
tion to the binding energy for all momenta. 
The nucleon momentum distribution 



n{p) 



^A{p,Lu)f{uj) 



(8) 



is very different different from the Fermi-Dirac distribu- 
tion. Momentum states below the Fermi momentum are 
depleted and a tail in the distribution n(p) for large mo- 
menta appears. The T-matrix approximation leads to a 
Fermi liquid behavior in the normal phase, with a jump in 



(6) the fermion density of Zp^ = ( 1 - 



.74 at the Fermi momentum. In the calculation the Fermi 
momentum is fixed by the constraint (|^) on the total den- 
sity. For a conserving approximation the Fermi momen- 
tum should be the same as the Fermi momentum of a free 
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Fig. 2. The binding energy for the T- matrix (solid line) and 
for the BHF (dotted line) calculations, and the Fermi energy 
for the T-matrix (dashed-dotted line) and for the BHF (dashed 
line) calculations as function of the Fermi momentum. 



fermion gas • This thermodynamical consistency re- 

lation is verified to a good accuracy by our calculations for 
the range of densities studied. 

The energy per particle, in the case of only two-body 
interactions, can be obtained from the single-particle spec- 
tral function 



Fig. 3. The quasiparticle energy uip and the average energy 
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+ u]A{p,u)S{uj) . (9) 



The binding energy per nuclcon as function of the Fermi 
momentum is presented in Fig. ^ for the self-consistent 
T-matrix approximation and compared to results from G- 
matrix calculations using the continuous choice of the aux- 
iliary potential. The results of the T-matrix approach lie 
above to the BHF binding energy for densities close to the 
phenomenological saturation point. Since we know that 
further hole line corrections do not modify the continuous 
BHF results drastically, we get an assessment of the accu- 
racy of the T-matrix approach. The higher the density the 
larger the discrepancy becomes. Correspondingly the sat- 
uration point in the T-matrix approach is shifted to lower 
densities (p = l-4po instead of 2.4po) and lower binding 
energies (the binding energy is reduced by 4MeV at po)- 
Very similar results are found for the equation of state of 
pure neutron matter . We note the the Hugenholz-Van 
Hove condition (^) is very well satisfied. 

The origin of the of the binding energy in the T-matrix 
approximation can be understood writing Eq. (ph as 



E/N IP 

' (271 



1 



2 \ 2m 



(10) 



LUp (hlh as function of momentum. 



with 



duj 
2^ 



i^MP,^)f{^)/n{p) , 



whereas in quasiparticle approaches it is 



E/N 



d^p 1 
- — r^nlp) — 
(27r)3 ^^'2 



P 

2m 



(11) 



(12) 



In Fig. 1^ we compare the removal energy to the quasi- 
particle energy LUp. Due to a large contribution of the 
background strength of the spectral function lying below 
the quasiparticle peak the removal energy is much below 
the quasiparticle energy. On the other hand, the single- 
particle energy tOp in the T-matrix approximation is gen- 
erally above the one obtained from G-matrix calculations. 
It can be seen by comparing the Fermi energies in the two 
approximations in Fig. ^j. These differences explain why 
the Fermi energy in the conserving T-matrix approxima- 
tion is equal to the binding energy at the local saturation 
point, following the Hugenholz-Van Hove relation (p. The 
binding energy is determined by uJ^ and Ep = LUp^ for the 
T-matrix calculation, whereas both quantities are deter- 
mined by LUp in the G-matrix scheme. 

This paper presents the first results self-consistent T- 
matrix calculation of saturation properties of nuclear mat- 
ter with a realistic potential. The binding energy obtained 
is smaller than the BHF result with the continuous auxil- 
iary potential. We note that a very similar shift in binding 
energy is observed in a BHF calculation, when including 
the rearrangement terms contribution to the binding en- 
ergy [ p^ . In Ref. | p5| by considering rearrangement term 
corrections to the single particle energies and to the bind- 
ing energy an improvement of the fulfillment Hugenholz- 
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Van Hove relation is found. The same can be observed in 
the T-matrix approach, we destroy a bit the binding en- 
ergies and improve considerably the single particle ener- 
gies from the BHF approach to get the relation (Q) right. 
We expect that after inclusion of ring diagrams contri- 
butions, as well as higher partial waves and three body 
forces corrections, the results on the saturation properties 
of nuclear matter of modern BHF approaches will be re- 
covered. The calculation of these corrections is standard 
and not related to the T-matrix approach. On the other 
hand, the real advantage of the self-consistent T-matrix 
approximation shows itself in the single-particle proper- 
ties. The quasiparticle energies lead to a Fermi energy 
consistent with the Hugenholz-Van Hove relation. The 
spectral function for large negative and positive energies 
can be calculated, with implications for the binding en- 
ergy and applications for electron scattering. We find at 
normal nuclear density an effective mass m* = .9m and 
Zp^ = .74. These numbers can be used in the estimation 
of reduced in medium cross sections and density of states 
at the Fermi energy. Generally, Zp does not change much 
with density or with the details of the potential used, pro- 
vided the lowest partial waves are described reasonably. 
We confirm the thermodynamical consistency of the nu- 
merical solution of the T-matrix scheme |^ for realistic 
interaction with several partial waves also. Finally let us 
note that the self-consistent T-matrix calculation can be 
straightforwardly extended to finite temperatures. 



16. P. Boiek, aucl-th/0202045 



17. P. Boiek and P. Czerski, |iiucl-th/020401S . 

18. J. Haidenbauer and W. Plessas, Phys. Rev C30, 1822 
(1984); ibtd, C32, 1424 (1985). 

19. J.M. Luttinger and J.C. Ward, Phys. Rev. 118 1417 
(1960); J.M. Luttinger, Phys. Rev. 119 1153 (1960). 



This work was partly supported by the KBN under 
Grant No. 2P03B02019. 



References 

1. K.A. Brueckner and J.L. Gammel, Phys. Rev. 109, 1023 
(1958); B.D. Day, Rev. Mod. Phys. 39 719 (1967) J. P. 
Jeukenne, A. Leugeunne and C. Mahaux, Phys. Rep. 25, 
83 (1976). 

2. R. Brockrnann and R. Machleit, Rhys. Rev. C42, 1965 
(1990). 

3. R.B. Wiringa, R.A. Smith and T.L. Ainsworth, Phys. Rev. 
C29, 1207 (1984). 

4. H.Q. Song, M. Baldo, G. Giansiracusa and U. Lombardo, 
Phys. Rev. Lett. 81 1584 (1998). 

5. W.H. Dickhoff, Phys. Rev. C58 2807 (1998); W.H. Dick- 
hoff, C.C. Gearhart, E.P. Roth, A. Polls and A. Ramos, 
Phys. Rev. C60 4319 (1999). 

6. P. Bozek, Phys. Rev. C59 2619 (1999). 

7. P. Boiek, Nucl. Phys. A657 187 (1999); P. Bozek, Phys. 
Rev. C62, 054316 (2000). 

8. Y. Dewulf, D. Van Neck and M. Waroquier, Phys. Lett. 
B510, 89 (2001). 

9. P. Bozek and P. Czerski, E ur. Phys. J. All, 271 (2001). 
10. P. Bozek, |nucl-th/0201084 



11. G. Baym, Phys. Rev. 127 1392 (1962). 

12. N.M. Hugenholz and L. Van Hove, Physica 24 363 (1958) 

13. M. Baldo, F. Bombaci, G. Giansiracusa, U. Lombardo, C 
Mahaux and R. Sartor,, Phys. Rev. C41 1748 (1990). 

14. F. de Jong and R. Malfliet, Phys. Rev. C44 998 (1991) 



15. P. Czerski, A. De Pace and A. Mohnari, |nucl-th/0204008 



